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Abstract 

Exact results stemming directly from Einstein equations imply that inhomogeneous Uni- 
verses endowed with vanishing pressure density can only decelerate, unless the energy den- 
sity of the Universe becomes negative. Recent proposals seem to argue that inhomogeneous 
(but isotropic) space-times, filled only with incoherent matter, may turn into accelerated 
Universes for sufficiently late times. To scrutinize these scenarios, fully inhomogeneous Ein- 
stein equations are discussed in the synchronous system. In a dust-dominated Universe, the 
inhomogeneous generalization of the deceleration parameter is always positive semi-definite 
implying that no acceleration takes place. 
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The homogeneity and isotropy of the geometry are not essential ingredients to establish 
a number of relevant results in relativistic cosmology. For instance, from the early sixties to 
the early seventies (see [1, 2, 3] and references therein), a research program on the singular- 
ity properties of general cosmological solutions has been conducted without relying on the 
isotropy and on the homogeneity of the geometry. 

The theme of the present paper is somehow opposite to the one analysed in [1, 2, 3] where 
the emphasis was on the role of the inhomogeneities (and anisotropies) in the proximity of a 
cosmological singularity. In the present paper we would like to understand if an inhomoge- 
neous space-time, filled with incoherent matter, can be turned into an accelerating Universe 
at later times. The inhomogeneities considered in the present investigation may arise dur- 
ing an early inflationary stage when quantum mechanical fluctuations of the geometry and 
of the inflaton field are inside the Hubble radius. Depending upon the parameters of the 
inflationary phase, the initial quantum fluctuations will be amplified leading to a quasi-flat 
spectrum of curvature perturbations that accounts, through the Sachs- Wolfe effect, for the 
tiny temperature ripples detected in the microwave sky by several experiments. While the 
temperature of the Cosmic Microwave Background (CMB) is in the range of few K degrees, 
the inhomogeneities are in the range of the /zK. This hierarchy of scales also implies that 
curvature fluctuations, say, after equality and before decoupling, are rather well described 
in the framework of linear theory. 

In recent months, it has been claimed [4] that since spatial inhomogeneities may have a 
non-trivial time dependence on length-scales larger than the Hubble radius, then, depending 
on the specific properties of the inhomogeneity, the deceleration parameter may well be 
negative, implying an effective acceleration of the Universe. While this proposal has been 
already criticized on various grounds [5, 6, 7, 8, 9] we find appropriate to scrutinize such a 
statement in the light of a set of fully inhomogeneous Einstein equations. Consider, indeed, 
Einstein equations 

t = t;, (i) 



R v „ = 8%G 

fj, 
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where R v is the Ricci tensor. In the specific case of a perfect relativistic fluid 

T» = (p + p)uX ~ P5 U », T; = (p-3p), (2) 
where is the velocity field obeying 

g^u,u u = 1, (3) 

while g^ u is the inverse of the four-dimensional metric tensor. The form of given in Eq. 
(2) (together with the synchronous condition) excludes the presence of torque forces. 

Being a symmetric rank-two tensor in 4 dimensions, it contains 10 independent entries, 
but 4 of them are associated with the freedom of choosing a coordinate system. Without 
loss of generality, the synchronous form of the line element may be adopted 

ds 2 = g^(x, t)dx p 'dx u = dt 2 — 7y(x, t)dx % dx 1 , (4) 
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where the symmetric 3-dimensional tensor / y i j(x,t) appearing in Eq. (4) carries 6 degrees 
of freedom which also correspond to the correct number of initial conditions for the general 
treatment of the problem. Equation (4) allows the determination of the components of the 
extrinsic and intrinsic (i.e. spatial) curvatures, 



K — K\, TrK 2 



K 3 K l 



Q -pm _ Q pm _|_ pmpf _ p£ pra 

' ij — U m L UjL mj tl ij 1 m £ 1 j m L i£ , 



(5) 
(6) 



where the three-dimensional Christoffel symbols are constructed using directly 7^. The (00), 
(ij) and (iO) components of Eq. (1) become, respectively: 



K - TrK 2 
1 d 



8ttG 



(p + p)u u c 



p- p 



7<9t 



yff K? ) - r? = &KG 



-(p + p)uiU J + ^—^Sl 



ViK - V k Kl = 8nG(p + p) Ui u°, 



(7) 
(8) 
(9) 



where r\ = 7^ r ki ; the overdot denotes a partial derivation with respect to t while V, denotes 
the covariant derivative defined in terms of 7^; in Eq. (8) 7 = det7ij satisfying 2K = —7/7. 
Equations (7), (8) and (9) have been used, in a slightly different form, in [1, 2, 3] to scrutinize 
the properties of inhomogeneous (and anisotropic) Universes. 
Using Eq. (4), Eq. (3) becomes 

u u° = 1 + Y j UiUj > 1, (10) 

where the second inequality holds since 7^ > 0. By definition of K and TrK 2 it is clear that 

1 



TrK 2 > -K 2 > 0, 
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where the first equality is reached in the isotropic case. Finally, the general definition of the 
deceleration parameter can be easily read-off from Eq. (7) and it is 



q(x, t) = 



K - TrfT 2 
TrK 2 



-1 + 



K 



TrK 2 



(12) 



In the isotropic limit j(x,t) — > a 2 (t)5ij, K\ — > — (d/a)5| and, from Eq. (12) q(x,t) — > 
—da/a 2 . Given the results reported in Eqs. (10) and (11), Eq. (7) implies that, since 
UqU° > 1, 



K - TrK' 



8nG 



(p + p)u u° + ^ 



> 47rG(3p + p), 



(13) 



implying, together with Eq. (11), that in the fully inhomogeneous case q(x,t) > provided 
(3p + p) > 0. If (3p + p) < the matter content of the Universe violates the strong energy 
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condition. No assumption has been made on the geometry as being split into a dominant 
background supplemented by a (comparatively small) inhomogeneous perturbation. 

Assuming now (as done in [4]) that the matter content of the Universe is provided by 
a dusty fluid with p — 0, the Universe can accelerate (i.e. q(x,t) < 0) only if the energy 
density p becomes negative ( violating simultaneously the weak, strong and dominant energy 
consitions). By taking the difference between Eq. (7) and the trace of Eq. (8) the standard 
form of the Hamiltonian constraint can be obtained, in the case p — 0, 

K 2 -TrK 2 + r = lQ7vGpu u°, (14) 

i.e. the negativity of the energy density also implies, according to Eq. (14), that K 2 — 
TrK 2 + r < 0. 

Before proceeding further, it is appropriate to note that Eqs. (7), (8) and (9) are sufficient 
to fully determine the dynamical evolution. However, it is always wise to keep an eye on 
the exact form of the covariant conservation equation whose two components become, in the 
case p = 

Id 1 

iVlpu u ) + —di(^f pu u l ) - pK ik u k u l = 0, (15) 



' d -[^pu°u i } + -± = d k ( y fy f m k u i ) -2pK l k u°u k + pT> fc u e = 0. (16) 



Every arbitrarily complicated solution of Eqs. (7), (8) and (9) must also satisfy, for consis- 
tency, Eqs. (15) and (16). 

We are now in condition of discussing some possible solutions of our system with the aim 
of showing that spatial inhomogeneities may turn the deceleration parameter from positive 
to negative. To do so the spatial inhomogeneities should be rather strong, so it will not be 
appropriate to separate the geometry into a homogeneous component supplemented by some 
small fluctuation. Hence, following Ref. [4] we can assume 

ltJ (x,t)=a 2 (t)e- 2 ^5 tJ , (17) 

that implies, according to Eq. (5), 

K\ = -(H -V)8{, K=-3(H-V), TrK 2 = 3(H - ^) 2 . (18) 

From Eq. (6), r\ and r are instead: 

r\ = — [V 2 ^ + did>V - (W) 2 5f + diV&iV], (19) 
a 2 

r = — [4V 2 ^-2(W) 2 ], (20) 

or 
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where V 2 ^ = 5 ij and (W) 2 = 5 ij d^d^. Equations (7), (9) and (14) lead then to 

the following system: 



3# - 3^ 2 + QH^ = AttG 



,2* 



p + 2(p + p)—u< 



3 2* 



d^f = AnG(p + p)Ui\ 1 + — u 2 , 



,2* 



^2 _ 2H ^ + [2V 2 ^ - (W) 2 ] 
3a 2 



8ttG 



P + (P + P)- 



,2* 



-U 



(21) 
(22) 
(23) 



where u 2 = 5^ Ui Uj and where the energy density has been separated, for later convenience, 
as p — p + p with p satisfying the usual Friedmann equations, i.e. 



H 2 = 



8ttG 



H + H = 



4ttG 



(24) 



and implying 2H = —3H 2 . Eliminating one of the p from Eq. (21) through Eq. (14) the 
equation to be solved becomes 



2^ - 3^ 2 + QH^f 



,2* 



8ttG 



3a 21 ' ' v ' ~ 7 J 3 
We can now look for a solution of the system in the form 

V(x,t) = f + \(t)e 2 fV 2 f, 



[2V 2 M>- (W) 2 ] = ^{p + pf—u 2 . 



2* 



(25) 



(26) 



where / = f(x) encodes the information on the large-scale inhomogeneities of inflationary 
origin and where e 2 ^V 2 / is the expansion parameter . In this expansion terms like e 2 -^(V/) 2 
will be considered subleading in comparison with e 2 ^V 2 /, i.e. |V 2 /| 3> (V/) 2 . Within this 
ansatz the momentum constraint can be solved, to lowest order as 



Ui = 



4irG(p + p)' 



(27) 



Plugging this solution back into Eq. (25) we find that the term containing u 2 is of higher 
order. The term \1/ 2 is also of higher order since it contains A 2 e 4 ^(V 2 /) 2 . The remaining 
terms cancel provided 



(28) 



If we now normalize the scale factor as a(t) = (t/t ) 2 ^ 3 (where t is the present time) A(t) is 
easily determined to be 



Ht) = ^-t 2 a(t) = 



T a(t), 



W'"'"' lhHf^ h (29) 

where, according to Eq. (24), we used Hq = 2/(3 to)- The solution given in Eqs. (26) and 
(29) coincides with the one given in Ref. [4]. Since terms (V/) 2 have been neglected in 
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comparison with V 2 /, a fortiori, in Eq. (25) the term 3\I/ 2 is subleading with respect to 2\l/. 
In fact \i/ 2 is of order A 2 e 4/ (V 2 /) 2 while * is of order \e 2f V 2 f. 
The obtained solution can now be parametrized as 

V = f(x) + a(t)Mx). (30) 

where the time-independent quantity \l/o is effectively treated, by the authors of Ref . [4] , as 
a free parameter. Inserting Eq. (30) into Eq. (12) we obtain 

3/2-a(t)fr /2 

g(M) = " 1+ (l-a(t)^o) 2 ' (31) 

that coincides with the expression of the deceleration parameter given in [4]. 

From Eq. (31), the authors of Ref. [4] conclude that : (i) if \&o — 0, then q(x,t) = 1/2; 
(%%) if |$o| < 1 ( an d < 0) q(x,t) < 0; (Hi) for large a(t) (i.e. t — > +oo) q(x,t) — > — 1. 
The only correct statement is (i). The other two are incompatible with the approximations 
made in solving both the dynamical equations and the constraints. 

In simple terms, the reason is the following. The deceleration parameter given in Eq. 
(31) is derived from Eq. (12) using the fact that TrK 2 = 3(H - ^) 2 = 3H 2 (l - aU ) 2 : 
but this is not correct since terms \I/ 2 have been neglected in the solution (because of higher 
order). Thus, if we treat as a free parameter we have to make sure, in particular, that 
from Eq. (25) the term 2\l/ is always leading with respect to 3^ 2 : using the parametrization 
(30) and the relation 2H = —3H 2 the requirement 2|\f r | > 3|\1/| 2 implies 

Ht)%\ < ^ (32) 

that is the most restrictive condition. If we now take into account that * = |e 2 -^V 2 /| < 1 
and that < 0, Eq. (32) implies an upper limit on a(t), i.e. 

a(t) < (33) 

implying that for finite is always bounded from above so that, for the smallest value of 
\l/o allowed in [4], i.e. ~ — 1/4, a < |. Therefore, the expansion only holds locally in 
time around t ~ to- For t ^> to higher gradients (and higher powers of the scale factor) 
must appear. This point has also been correctly emphasized in [7] and we checked it both 
by going to higher order in perturbation theory and by extending the present analysis to the 
other degrees of freedom present in 7^ together with \P. 

In the limit a — > 00 the terms neglected in the solution given in Eqs. (26) and (29) 
become important and the momentum constraint is not satisfied anymore. The limit a — > 00 
can be only taken, according to Eq. (33), if ^0 ~~ 0- But, in this case q — 1/2. 

To complete the previous series of statements the correct evaluation of the decelera- 
tion parameter will now be discussed. Since TyK 2 (appearing in the denominator of Eq. 
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(12)) must be positive semi-definite, its expression, within the approximations made in the 
derivation of Eq. (30), is 

TrK 2 = H 2 (l-2cM ), (34) 

with ty^a < 1/2; the latter condition is less restrictive than the one of Eq. (32). Therefore, 
using the observations reported above (and in particular of Eq. (34)), the deceleration 
parameter can be determined from Eq. (12) as 

Taking now into account Eq. (32) we can easily see, from Eq. (35), that q(x, t) as a function 
of a^o is always positive semi-definite. It is amusing that q(x,t) — > for the the largest 
(negative) value of ^ a, i.e. —1/3. This strongly suggests that the large time limit of q(x,t) 
for a maximally underdense Universe (still compatible, though, with perturbation theory) is 
as argued in a specific (but exact) example by the authors of Ref. [7]. 
Since the energy density is, to leading order in |e 2 ^V 2 /| 

P = P + P = ^[l + Sa%], (36) 

when q(x,t) — > 0, p and p exactly cancel (as it follows, in the case p — 0, from Eq. (13)). 
Thus, for |a\& | < 1/3, p > and all the energy conditions are correctly enforced. Equation 
(36) may also be derived from Eq. (15). To leading order in |e 2 ^V 2 /| (neglecting terms 
(V/) 2 and implementing the momentum constraint), Eq. (15) implies: 

■p + 3Hp = 0, 'p + 3Hp - 3^p = 0. (37) 

Equation (36) is solution of Eqs. (37) provided terms going as \I/ 2 are neglected. If this is 
not the case covariant conservation equations and Einstein equations do not lead any longer 
to compatible solutions. Once more, if the limit a — > oo is taken blindly in Eq. (36), p 
becomes arbitrarily negative when ^ < (underdense regions). But this procedure is not 
consistent since, at most, \a^o\ < 1/3 and |^ | < 1- 

In [5] the authors argue, correctly, that a renormalization of the local spatial curvature 
cannot imply that the Universe is accelerating. Our result shows that this is indeed the 
case: the acceleration cannot be obtained unless the inhomogeneous solution is extrapolated 
in a regime where the terms neglected in the perturbative expansion become dominant. As 
previously remarked, also the results of [7] are consistent both with the present treatment 
and with the criticism raised in [5]. The authors of Ref. [7] select a specific f(x) that 
leads to an exactly solvable model with asymptotically vanishing acceleration parameter. 
By comparing the exact result with the proposal of Ref. [4] , the authors of Ref. [7] are led 
to conclude that the apparent acceleration is just a consequence of neglecting higher orders 
in the gradient expansion. This conclusion fits well with our findings. 
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Criticisms of the proposal [4] have been also presented in [6, 8, 9]. While these papers 
raise important issues, they are more related with observational implications of a class of 
solutions whose debatable correctness was the aim of the calculations presented here. The 
interesting points raised in [6, 8, 9] together with generalizations of the arguments presented 
here will be the subject of a forthcoming publication. 

To conclude, the main lessons to be drawn from the present investigation are the following: 
(a) in the context of the perturbative expansion proposed in [4] the deceleration parameter of 
a matter dominated Universe is always positive to a given order in perturbation theory; (b) 
the energy density is always positive semi-definite to a given order in perturbation theory; 
(c) the claimed acceleration is the result of the extrapolation of a specific solution in a regime 
where both, the perturbative expansion breaks down and the constraints are violated. 

References 

[1] E. M. Lifshitz and I. M. Khalatnikov, Sov. Phys. Usp. 6, 495 (1964). 

[2] V. A. Belinskii, E. M. Lifshitz, and I. M. Khalatnikov, Sov. Phys. Usp. 13, 745 (1971). 

[3] I. M. Khalatnikov and E. M. Lifshitz, Phys. Rev. Lett. 24, 76 (1970). 

[4] E. Kolb, S. Matarrese, A. Notari and A. Riotto, arXiv:hep-th/0503117. 

[5] G. Geshnizjani, D. J. H. Chung and N. Afshordi, arXiv:astro-ph/0503553. 

[6] E. Flanagan, arXiv:hep-th/0503202. 

[7] C. M. Hirata and U. Seljak, arXiv:astro-ph/0503582. 

[8] S. Rasanen, arXiv:astro-ph/0504005. 

[9] E. Siegel and J. Fry, arXiv:astro-ph/0504421. 



7 



